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Abstract: It is known that the entanglement entropy of a scalar field, found by tracing over 
its degrees of freedom inside a sphere of radius 1Z, is proportional to the area of the sphere 
(and not its volume). This suggests that the origin of black hole entropy, also proportional 
to its horizon area, may lie in the entanglement between the degrees of freedom inside and 
outside the horizon. We examine this proposal carefully by including excited states, to check 
probable deviations from the area law. 

PACS Nos.:04.60.-m,04.62.,04.70.-s,03.65.Ud 

Resume : Nous savons que l'entropie d'entrelacement d'un champ scalaire, trouvee en suivant 
ses degres de liberte a l'interieur d'une sphere de rayon 1Z, est proportionelle a la surface 
de la sphere (et non a son volume). Ceci suggere que l'origine de l'entropie d'un trou noir, 
egalement proportionelle a la surface de son horizon, peut se trouver dans l'entrelacement 
entre les degres de liberte a l'interior et a l'exterior de l'horizon. Nous examinons avec soin 
cette hypothese, en incluant les etats excites, afin d' identifier les deviations possibles de la loi 
de surface. 

[Traduit par la redaction] 



1. Introduction 

There are strong indications that a black hole (BH) of mass M and horizon area Ah possesses 
entropy and temperature, given respectively by [ 1 ] : 

Sbh = ^, Th =8^GM' Planck length). (1) 

The above entropy and temperature satisfy the laws of BH Thermodynamics 

Th = Constant on horizon , d (Mc 2 ) = TndSBH + work terms . A(Sbh + S ma tter) > (2) 

where 'work terms' are relevant for BHs with charge(s) and angular momenta, and S ma tter refers to 
the entropy of matter outside the BH. Usually the latter can be given a microscopic interpretation by the 
relation S ma tter = l n ^> where is the number of accessible micro-states, compatible with the various 
macroscopic parameters such as temperature, pressure and volume. Although such an interpretation for 
Sbh is incomplete, important progress has been made in various approaches to quantum gravity [2-5]. 
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This, coupled with the area proportionality of BH entropy (as opposed to volume proportionality), has 
raised a fundamental question in quantum gravity, namely: What is the origin ofBH entropy? Potential 
candidates include strings, D-branes, spin-network states and conformal degrees of freedom at the 
horizon. Another popular (and also incomplete) approach to the BH entropy is entanglement entropy, 
which is a measure of the information loss due to the spatial separation between the degrees of freedom 
inside and outside the horizon. The so-called brick wall model has been a concrete realization of this 
idea, in which entanglement entropy of the scalar (and other) fields have been computed by tracing over 
their degrees of freedom outside the horizon [6] 1 . The brick-wall entropy turns out to be proportional 
to the BH horizon area. The problem, however, is that due to the infinite growth of the density of 
states close to the horizon, one has to impose ultraviolet cut-off near the horizon and hence the entropy 
depends on the cut-off scale. This clearly is an undesirable feature. 

A simpler physical system was considered in refs. [7, 8], in which the entropy of a scalar field on 
a suitably discretized flat space was computed numerically, by tracing over the degrees of freedom 
of a hypothetical sphere of radius 1Z. This gave the remarkable result that the entanglement entropy 
was indeed proportional to the area of the spherical surface (also see [9] for analytical proofs of the 
result). This further supported the idea that entanglement was responsible for BH entropy. One of 
the key assumptions in refs. [7, 8] was that the harmonic oscillators (HOs) resulting from the scalar 
field, on discretization of space, were all in their ground states. In this article, we would like to relax 
this assumption, and investigate the robustness of the entropy-area relation in these physical systems. 
As a step towards understanding of the entanglement entropy of N coupled harmonic oscillators (by 
tracing over n < N degrees of freedom), in this work we consider two HOs (N = 2, n = 1) in two 
physically interesting limits: coherent states and superposition of excited and ground states. We show, 
numerically, that the presence of excited states results in an increase of entropy. The generalization to 
arbitrary N will be left to a future publication [10]. 

In the next section, we will briefly review the entanglement entropy for ground state HO. In section 
(3), we will generalize the results for the two HO wave-function which is a superposition of ground 
state and the first excited state, and show that the entropy increases. In the concluding section, we will 
remark on the possible implications of our results to the physically interesting case of a large number 
of oscillators, whose couplings are determined by the Lagrangian of a free scalar field. 

We will follow the notations of ref. [8] to provide easy comparison, and henceforth use Ti = 1. 



2. Ground State Entanglement Entropy 

The Hamiltonian for two coupled HOs of unit masses: 
1 



H =2 



P1+V2 + ko (x 1 + x 2 ) + h (x 1 - x 2 y 



(3) 



has the ground state solution in terms of normal modes (x± — Xl ^ 2 , ^+ — Vh), U- = V^o + 2fci): 



tpo (xi,x 2 ) = ii (x+)ipo{x-) 



1/4 



r l/2 



■ exp [— (ijJ+x 2 + + oj-x 2 J) /2] . 



When traced over the oscillator characterized by x\, the resultant density matrix is: 

/oo 
dxiipa (x 1 ,x 2 )ipo (xi,x' 2 ) = 
-oo 



P0 



7-/3 



cxp [-7 (x 2 2 + x'i) /2 + (3x 2 x' 2 ] 



(4) 



(5) 



Note that tracing over the outside does not pose any conceptual problem, since for a pure system, tracing over a given 
subsystem and its complementary subsystem yield identical entropies [11]. 
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whose eigenfunctions and eigenvalues, respectively are 

f n (x) = H n (yftix) exp (~ax 2 /2) , p n = (1 - £) C • 

The above density matrix gives rise to the following entropy [7, 8]: 

OO £ 

S(0 = -Tr (p In p ) = - ]T Pn ln Pn = - ln(l - £) - 737 In £ , 



71=0 



where 



i? 2 = 



2\2 



< 1, a = uj-R, f3 = 



u;_(l-^ 2 ) 
4(1 + R 2 ) 



7 = 



1 + 6R 2 + R 4 
4(1 + i? 2 ) 



l-R 
1 + R 



(6) 



(7) 



(8) 



Note that R = 0, 1 correspond to the strongly coupled and uncoupled limits respectively. 
The Hamiltonian for a free, massless scalar field ip in flat space-time is given by: 



H=ljd 3 x [7r 2 (r) + |W)| 5 



(9) 



Discretizing the space (a being the lattice spacing, Na signifying the infrared cutoff, and I, m are the 
parameters in the spherical harmonics Y^ m (6, </>)): 



1 N \ ( 1\ 2 / 



<Plm,j_ __ Pirn J + 1 

i + i 



+ — ^—<Pi 



(10) 



Ira j — 1 _ 

The above expression for iJ; m is a special case of the general TV-coupled oscillator Hamiltonian: 

N N 



The corresponding iV-HO ground state wave function is given by: 



ipo(xi,. 



ini 



-,1/4 



exp 



a; ■ O • a; 



(11) 



(12) 



where il 2 = K. It can be shown that for the ground state wave function, the density matrix (when one 
traces over n < N oscillators) can be factorized into a product of (N — n) 2-HO density matrices. Thus 
the total entropy is simply the sum of the entropies. For the scalar field, n is taken to be proportional 
to the radius of the sphere which is traced over, i. e., 1Z = (n + l/2)a. For the Hamiltonian (10) the 
interaction matrix Kij can be read-off, resulting in the entanglement entropy 



S = 0.3(n+ 1/2) 2 <xK 2 , 



(13) 



signifying area proportionality. It is worth noting, again, that all the oscillators are assumed to be in 
their ground states. 

3. Entanglement entropy of excited states: 

Let us now consider the excited states of the N HOs discussed in the previous section. The corre- 
sponding wave-function is: 



ip(xi, . . . ,X n ,X n+ i, . . .,X N ) = 



1/4 



exp 



x ■ Ct ■ x 



i — l v 



(14) 
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where Kb = UKU T is a diagonal matrix (U T U = In), x = Ux and = 1 ... N) are indices of 
the Hermite polynomials. The density matrix, tracing over first n of N oscillators, is 



PO (x n 



-i-i, ... , xn', x n+ i, 



N 

n 



\9.\ 



-,1/2 



rJV 



Y[dx i exp 



• Q, ■ x 



H Vi K D x, exp 



■ n-x' 



N 1 

TT ff L (K D ix'). 



(15) 



The evaluation of the integral of the product of 2N Hermite polynomials, although may not be im- 
possible, is in general, non-trivial. In order to keep the calculations simple, we consider two specific 
physical cases: (i) coherent states and (ii) superposition of ground and first excited states. 

The coherent states, which are eigenstates of the harmonic oscillator annihilation operator with real 
eigenvalues, are described by the following wave function: 



ipcs(x, a)=^ {x-a) = e lpa ip (x) 



(16) 



The expectation of the position operator, w.r.t the coherent state wave function, oscillates in time with 
an amplitude a and the state has the minimum allowable uncertainty 



Ap Ax = - 
2 



(17) 



same as that of the ground state. For two coupled oscillators, the corresponding coherent state is: 

i>cs{xi,x 2 ) = ipcs(x+,a)tpcs{x-,b) = ip {x+ - a)ip {x- - b) . (18) 

Defining x 2 — X2 — (a — b) j\pl , it is easy to show that the corresponding density matrix retains the 
same form as (5), albeit in terms of these new variables: 



PCS(X2,X' 2 ) = p out (X 2 ,X 2 ) 



(19) 



Thus, from Eqs. (6), it follows that the eigenfunctions are f n (x) and eigenvalues remain unchanged 
ip n ), and we get the interesting result that the entropy is the same as that for the ground state! Presum- 
ably, this is because of the fact that coherent states are obtained by translating the ground state in phase 
space. The result can be easily generalized to an arbitrary number of harmonic oscillators [10]. 
Next, we consider the superposition of the ground and first excited state of the 2-HO system: 

ip(xi,x 2 ) = a x ■0i(x + )-0o(a;-) + Pi "4>o(x+)ipi(x-) + 71 tp (x+)tpo{x-) [a\ + p\ + = 1],(20) 

,1/4 l 



where ^ n (x) = N n H !/2 R n (V^ x), N n (w) - (^) 



(21) 



is the n th excited state of an oscillator. Although from the identity of particles one would expect 
a\ = Pi, we do not impose such a condition at this point. From (15), the density matrix follows: 

p{x2, x' 2 ) = Po(x 2 , x' 2 )[A (x 2 + x 2 ) + B x 2 x 2 + C (x 2 + x 2 ) + D] , (22) 

where po(x 2 , x' 2 ) is the ground state density matrix given by Eq. (5), and the constants are given as: 

A = a\a + P\a3 + a.\P\a± , B = a\b + Pfb 3 + a\P\b\ , 
D = ale + P\c z + a 1 p 1 c i + , 



a 6 



I? 



C = 7i(aia 6 + Pia 7 ) , 



1 



R 2 {l-R 2 ){Z + R 2 )< 
4(1 + R 2 ) 2 



R 2 (5 



2R 2 



a 7 



1 + R 2 



i? 4 )co_ 



R 2 



2(1 + i? 2 ) 2 



1 + i? 2 ' 



(23) 
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(1- R 2 )(1 + 3R 2 )uj- 



as = - 



4(1 + R 2 ) 2 
1 - R 2 \ 2 lu_R 



&3 = 



(1 + 2R 2 + bR 4 )u- 
2(1 + R 2 ) 2 



1 



C3 



1 + R 2 J 2 ' 



&4 = - 



R(l + 6R 2 + R 4 )u 
(1 + i? 2 ) 2 



1 + R 2 ' 
2R 



C 4 = 



1 + i? 2 



It can be verified that: Tr(/>) = f_*^ cZx 2 p(x 2 ,x 2 ) = a 2 + {3 2 + -f 2 = 1 . To find the eigenvalues of the 
density matrix (22), we follow the general procedure outlined in ref. [12]. First, we expand p(x 2 , x 2 ) in 
terms of general HO eigenstates (although, in principle any complete set of functions should suffice): 

p(x 2 ,x 2 ) = 22 h n (x 2 )g n (x 2 ) , h m (x 2 ) = N m (a) exp ( — J H m (y/a x 2 ) . (24) 

Inverting, we get: 



i(x 2 ) = / dx 2 p(x2, x 2 ) h m (x 2 ) (25) 

J — oc 

= p m N m e [(B 1 ^ + £; 1 )^ m+1 (^^) + (c , ia;' 2 2 + I) 1 +Fia;' 2 )/f m (^a;2)] , 



where 

Si = -Va 



207 



+ 



7 2 - a 2 ^ 7 2 _ Q 2 



+ c, D12 



2aa 
7 Z — a z 



7 — a y 7 — a 



£>i = Dn +Di 



F 1= d 



1 + 



7 + a 



7 — a 



(26) 



7 + a 7^ - a^ ' ^2 - a 2 

a = a\a + p\az + ai/3i04, 6 = a 2 b + /3 2 63 + ai/?i&4, c = a 2 c + /3 2 C3 + ai/?iC4 + 7 2 , d = aide + /?ifl7 , 
and dj, 6, have been defined in Eq. (23). The next step is to define the matrix equivalent of p, i. e., 



/oc 
dx g„ 
-OO 



(x)hp(x) 



(27) 



D n +pD 12 + 



Bi(p+1) , Ci(2p + 1) 



2a 



^ 

S pm + -i\/(p+l)(p + 2)5 p ,, 



2a 



+ - 1) ( 4= + ^ J <W+2 + F iJ^- S p+i,™ + (Si V2p + Fi^^Vi.m) 5„-i 



\/a 

Although formally diagonalizable, the eigenvalues A p of the above penta-diagonal matrix are most 
easily found numerically. With MAPLE, using upto 40 x 40 matrices, we verified that it has unit trace. 



m=0 



Bi C\\ / B 1 Ci 

£i + ^ + - L +mL>i 2 + ^i + — 



The corresponding entropy as function of a\ , (3i , R defined as: 

00 

S(ai,(3i,R) = - ^ \ p In A p 

p=0 



•(28) 



(29) 



was also computed numerically, and for all ot\,(3\ ^ it was found that S(a\,(3i,R) > 5(0, 0, i?), 
where 5(0, 0, R) is the ground state entropy. The equality holds only in the uncoupled limit R = 1 and 
ai = These features are visible in Fig.(l), where we have plotted entropies for the excited state 



NRC Canada 



6 



Can. J. Phys. Vol. 




0.2 0.4 0.6 0.8 

R 



Fig. 1. Plots of the entanglement entropy of the excited state 5(l/\/2, R) (black curve) and that of the 
ground state 5(0, 0, R) (grey curve) vs. R. Note that 5(1/V2, l/y/2, R) > 5(0, 0, R) for all R < 1. 



[ai = Pi = 1/ \/2, 7i = 0] as well as the ground state. In brief, any amount of excited state in the su- 
perposition increases the entropy. This is intuitively expected, since it can be shown that the expectation 
of energy in the state (20) is give by: < if (ai, /3i, 71) >= ^ [a\(l + 3R 2 ) + /3 2 (3 + R 2 ) + 7?(1 + R 2 )] , 
from which it follows that < H{a 1 ,(3 1 ,^ 1 ) > - < ff(0,0,l) >= ^ {a x R 2 + (3 2 ) > . That is, 
the expectation of energy is least for the ground state; and higher energies are normally associated with 
higher entropies. 

4. Summary and Outlook 

In this paper, we have shown that the entanglement entropy of two coupled HOs, with coordinates 
of one traced over, is more for excited states compared to when they are both in their ground states. 
We would like to extend our results to N oscillators, with n < N of them being traced over. This 
would enable one to compute the entanglement entropy of a free scalar field in flat space-time when 
its degrees of freedom inside a given region are traced over, and check whether it is proportional to the 
area of the bounding surface [10]. It would also be interesting to extend the results to BH space-times 
with the surface mentioned above coinciding with its event horizon. We hope to report on these and 
related issues in the near future. 

Note added: After the submission of this paper to the Canadian Journal of Physics, we extended 
the work (reported here) to the free scalar field in flat space-time in Ref. [10]. We have shown that the 
entanglement entropy is proportional to area when the scalar field degrees of freedom are in generic 
coherent states, and first excited state, although in the latter case, the entropy increases manyfold. 
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